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ABSTRACT 

We present a detailed analysis of Chandra X-ray observations of the lensing cluster of galaxies 
CL0024+17 at z = 0.395. We found that the radial temperature profile is consistent with being isother- 
mal out to ~ 600 kpc and that the average X-ray temperature is 4.47i !s4 keV. The X-ray surface 
brightness profile is represented by the sum of extended emission centered at the central bright elliptical 
galaxy with a small core of 50 kpc and more extended emission which can be well described by a spherical 
/3-model with a core radius of about 210 kpc. Assuming the X-ray emitting gas to be in hydrostatic 
equilibrium, we estimated the X-ray mass within the arc radius and found it is significantly smaller than 
the strong lensing mass by a factor of about 2-3. We detected a strong redshifted iron K line in the 



X-ray spectrum from the cluster for the first time and find the metal abundance to be 0.76 



+0.37 
-0.31 



iolar. 



Subject headings: galaxies: clusters: individual (CL0024- 
lensing — dark matter 



17) — X-rays: galaxies — gravitational 



1. INTRODUCTION 

CL0024+17 is one of the most extensively studied dis- 
tant clusters of galaxies, located at z = 0.395. The cluster 
is known to exhibit spectacular multiple arc images of a 
background galaxy produced by the gravitational lens ef- 
fect. This makes the cluster a unique target for studying 
the distribution of not only luminous matter but also dark 
matter. 

The lensed arc system was initially noted by Koo (1988) 
and the first spectroscopic observation placed limits on the 
arc redshift, 1 < z arc < 2 (Mellier et al. 1991). From 
the deep HST imaging, four arc images of the blue source 
galaxy were identified around the tangential critical curve 
and further was a radial arc of the same source (Smail et 
al. 1996; Colley et al. 1996). In addition, CL0024+17 is 
the first cluster where a coherent weak shear signal was 
detected (Bonnet et al. 1994). 

The remarkably well-resolved lens images allow detailed 
modeling of the dark matter distribution in the cluster 
center and several authors have attempted the mass re- 
construction (Kassiola et al. 1992; Smail et al. 1997; Tyson 
et al. 1998; Broadhurst et al. 2000). Tyson et al. (1998) 
has built a very detailed mass map comprising 512 free 
parameters. They found that, excluding the mass associ- 
ated with individual galaxies, the dark mass distribution 
is smooth and has little asymmetry, and the central pro- 
file has a soft core of 35/i _1 kpc. This is contradictory to 
the central cusp predicted by the CDM numerical simula- 
tions (Navarro et al. 1996, hereafter NFW) and thus has 
motivated investigations of different types of dark matter 
such as self-interacting dark matter (Spergel & Steinhardt 
2000; Moore et al. 2000). Broadhurst et al. (2000) mea- 
sured the arc redshift to be z arc = 1.675 and also built 
a lens model in a simplified manner assuming NFW pro- 



files for cluster galaxies. Their model needed to include 
only the brightest eight cluster galaxies to reproduce the 
observed lensed configuration. They reached the conclu- 
sion that, in contrast to Tyson et al. (1998) 's model, the 
average mass profile is consistent with an NFW profile. 
Shapiro & Iliev (2000) pointed out that the fit by Broad- 
hurst et al. (2000) implies a cluster velocity dispersion that 
is much larger than the value measured for this cluster by 
Dressier et al. (1999). 

The galaxy velocity dispersion of a ~ 1200 km/s was 
found by Dressier & Gunn (1992) and Dressier et al. (1999) 
based on the data of about 100 galaxies. Czoske et al. 
(2001, 2002) constructed a galaxy catalog of the wide- 
field spectroscopic survey of the cluster, and newly identi- 
fied a foreground and a background group of galaxies well 
aligned along the line of sight. They clearly demonstrated 
the bimodal distribution of galaxies in the redshift his- 
togram for 300 objects in the neighborhood of CL0024+17. 
If restricted to the galaxies in the main component, the ve- 
locity dispersion is ~ 600 km/s. Thus if interpreted as a 
value for a relaxed cluster, the dynamical mass would have 
only a quarter of that previously derived (Schneider et al. 
1986). 

X-ray emitting intracluster gas is an excellent tracer of 
the dark matter potential. Soucail et al. (2000) performed 
a combined analysis of the ROSAT and the ASCA data 
and estimated the cluster mass within the arc radius un- 
der the assumption of hydrostatic equilibrium (the X-ray 
mass, hereafter). They found that there is about a factor 
of ~ 3 discrepancy between the X-ray mass and the strong 
lensing mass (Tyson et al. 1998; Broadhurst et al. 2000). 
They compared the extrapolated X-ray mass to the weak 
lensing mass within 3 h^ Mpc and found it is again lower 
by a factor of ~ 3. Because the ROSAT HRI image sug- 
gested an elongated gas distribution (see also Bohringer 
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et al. 2000), they considered that the discrepancy may be 
partly caused by the irregular mass distribution. 

The 'mass discrepancy problem' between X-ray and 
strong lensing mass estimation has been reported in many 
other lensing clusters and the X-ray mass is systemati- 
cally lower than the strong lensing mass by a factor of 2-5 
(Miralda-Escude & Babul 1995; Wu & Fang 1997; Ota et 
al. 1998; Hashimotodani 1999). A variety of possible ex- 
planations have been proposed; non-thermal pressure sup- 
port, effect of complex mass distribution, etc. (e.g. Hattori 
et al. 1999). However, there is as yet no definitive ex- 
planation, and it is possible that some important physical 
process may not have been considered in regard to the 
central region of clusters. 

For the specific case of CL0024+17, however, there were 
still large measurement uncertainties in both the X-ray 
temperature and the image morphology, due mainly to 
heavy contamination from a bright seyfert galaxy that lies 
near the cluster center. Thus for the cluster mass es- 
timation the temperature determination is crucial. The 
high-resolution Chandra data eliminate the contamination 
from point sources and enable us to determine the X-ray 
spectrum and the spatial distribution of the cluster gas 
simultaneously. In this paper, we report accurate mea- 
surements of the temperature and the morphology with 
Chandra, from which we consider the discrepancy between 
the X-ray and the strong lensing masses. 

We use H = 50 km/s/Mpc and fl = 1, and thus 1' 
corresponds to 383 h§ kpc at z — 0.395. The quoted er- 
rors are the 90% confidence range throughout the paper, 
except where noted. We use the solar abundance ratio of 
iron atoms to hydrogen atoms, Fe/H = 4.68 x 10~ 5 (An- 
ders & Grevesse 1989). 

2. OBSERVATION AND SOURCE DETECTION 
2.1. Chandra observation of CL0024+17 

We observed CL0024+17 with the Chandra Advanced 
CCD Imaging Spectrometer (ACIS-S) detector on Septem- 
ber 20, 2000. The target was offset from the ACIS-S nomi- 
nal aimpoint with a Y offset of — 1'.33 in order to avoid the 
outskirts of cluster diffuse emission falling in a gap between 
the CCD chips. The CCD temperature was -120°C. The 
data were reduced using CI AO version 2.2 with CALDB 
version 2.15. We removed periods of high background lev- 
els exceeding 3ct above the mean quiescent background 
rates. The net exposure time was 37121 sec (93.2% of the 
total exposure). To improve the Chandra astrometry we 
used the Aspect Calculator and corrected the aspect off- 
sets. In Figure 1, we show the ACIS-S3 image in the 0.5-7 
keV band. The strongest X-ray peak of CL0024+17 is at 
00:26:36.0, +17:09:45.9 (J2000) and extended emission is 
detected out to ~ 2' in radius. We refer to the X-ray peak 
as the Gl peak hereafter. 

We planned to set the roll angle at 127° with a tore- 
lance of ±8° in order to place the position of the "pertur- 
bation", indicated by the weak lensing analysis (Bonnet 
et al. 1994), on the ACIS-I2 chip. The perturbation is 
about 6'. 7 north- northeast of the gravitational shear co- 
herent to the CL0024+17 center. However the tolerance 
was exceeded by 8'. 4 and the actual roll angle was 135.14°, 
and thus the perturbation was outside the ACIS field of 
view. 



2.2. Source detection in the ACIS-S3 

We searched for point-like sources in the ACIS-S3 field 
with the WAVDETECT algorithm with a significance 
threshold parameter of 10~ 6 and detected 38 sources. In 
the following analysis these point sources were excluded 
with a radius of 7 times the size of the point spread func- 
tion (PSF) at the source position. The PSF size is defined 
as the 40% encircled energy radius at 1.5 keV. We cross- 
correlated the source positions with the ROSAT HRI and 
found 5 sources in the ACIS-S3 field of view (SI, S2, S3, 
S4, S5 in Soucail et al. 2000) were consistent with each 
other within 1" ~ 3", which is smaller than the point- 
ing accuracy of the ROSAT HRI (la error is typically 
6"; Briel et al. (1997)). For SI, there were two optical 
counterparts within the ROSAT error box (16" x 16"), 
which are a foreground starforming galaxy and one of the 
cluster member galaxies. The Chandra imaging clearly re- 
solved SI into two point sources, whose J2000 coordinates 
are respectively determined to be (a, <5)— (00:26:31.691, 
+ 17:10:21.71) and (00:26:31.101, +17:10:16.48). These are 
cataloged by Czoske et al. (2001) as #282 (z = 0.2132) 
and #267 (z = 0.4017) respectively, thus by comparing 
the positions we confirmed that the Chandra astrometry 
is accurate down to 0".3. 

2.3. MG0023+171 

The radio source MG0023+171 with z = 0.946 (He- 
witt et al. 1987) was covered with the ACIS-S0 chip dur- 
ing the observation. MG0023+171 has two optical coun- 
terparts separated by 5" which have been interpreted as 
gravitationally lensed images. The large separation an- 
gle implies a large mass-to-light ratio for the lensing mat- 
ter however the lensing object has not yet been identified. 
We searched for X-ray emission from a "dark object" , but 
did not detect any object using WAVDETECT with the 
threshold parameter of 10 -5 . The photon counts in the 
0.5-7 keV band within a circle of radius 5" centered at 
image A of MG0023+171 system is 6, without subtracting 
background. We thus constrained the 3cr upper limit on 
the X-ray energy flux as 8.3 x 10~ 15 erg/s/cm 2 assuming a 
power-law spectrum with the photon index of 2.0. On the 
other hand if we extrapolate the B magnitude of the radio 
source, 21.9, utilizing a typical aox index for radio-loud 
quasars of 1.6 (Green et al. 1995), the X-ray flux at 2 keV 
is expected to be 6.4 x 10~ 16 erg/s/cm 2 /keV. We found 
that this is below the detection limit of the current obser- 
vation. If we attribute the entire X-ray flux to a lensing 
object at an assumed redshift of 0.4, the upper limit of 
the X-ray luminosity is constrained to be ~ 7 x 10 42 erg/s. 
Thus a massive object such as a galaxy group or a galaxy 
cluster is unlikely to be the dark lens candidate as long as 
one assumes the normal X-ray properties. 

3. SPECTRAL ANALYSIS 

3.1. Overall spectrum 

We extracted the cluster spectrum from a circular region 
of r — 1'.5, centered on the Gl peak. The background 
was estimated from the 2'. 5 < r < 3'. 2 ring region and 
subtracted from the above spectrum. Note that we inves- 
tigated the positional dependency of the background by 
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comparing two spectra accumulated in the corresponding 
detector regions for the source and the background spec- 
tra of the blank-sky observation data. We found that the 
difference of the background spectra from the two regions 
is sufficiently small (the difference in intensity is less than 
20% in each spectral bin and the overall normalization de- 
viates by 4%) and the results of the present analysis are 
not affected within the statistical errors. We generated 
the telescope response file (i.e. the ARF file) for an ex- 
tended source with the WARF procedure, which sums the 
ARFs according to the weight of counts in each bin in the 
given image region. We found that the weighted ARF dif- 
fers only by 2% at most in the ACIS-S energy band from 
that made for a point source. Thus effect of telescope vi- 
gnetting is negligible in comparison to the statistical errors 
of the present data. 

We fitted the cluster spectrum in the energy range of 
0.5 — 7 kcV with the MEKAL thin-thermal plasma model 
(Mewe et al. 1985, 1986; Kaastra 1992; Licdahl ct al. 
1995). If we include the hydrogen column density of 
the neutral absorption as a free parameter, we obtained 
A H = (5.3 ± 2.5) x 10 20 cur 2 . This is consistent with the 
Galactic value, thus we fixed it at Ah = 4.2 x 10 20 cm~ 2 
(Dickey & Lockman 1990). The result of the fit is shown in 
Figure 2 and Table 1 . This provided a good fit to the data 
and the resulting \ 2 value was 47.5 for 45 degrees of free- 
dom. The temperature was determined to be kT = 4.47 
kcV with the 90% error range of 3.93 - 5.30 keV. The 
result is consistent with the previous measurement with 
ASCA, which had a large (~ 50%) uncertainty (Soucail et 
al. 2000). 

We detected strong line emission at about 4.8 keV 
in our reference frame, which is consistent with a red- 
shifted, highly ionized Fe-K line emitted from an object 
at z = 0.395. This is the first detection of the iron line 
from CL0024+17. We obtained the metal abundance to 
be 0.76toJi solar (90% error). In the Chandra spectrum 
there is also a significant contribution from the Fe-L com- 
plex below 1 keV. If we mask the energy bins between 4 
keV and 5.2 keV that cover the energy of Fe-K line, the 
MEKAL model fitting resulted in the metal abundance of 
1.38±J |f solar (90% error). The best-fit value is higher 
than that we obtained from the above fit however still 
consistent within the error ranges. Thus the metallicity 
for CL0024+17 is about a factor of 2 higher than those of 
typical distant clusters (Mushotzky & Loewenstein 1997; 
Matsumoto et al. 2000; Ota 2001). On the other hand, the 
previous ASCA spectral analysis by Soucail et al. (2000) 
did not constrain iron emission lines because the photon 
statistics of the cluster spectrum were seriously limited 
due to the contamination from the SI emission. The high 
metallicity value we derived is discussed later in the paper 
(See §6.1). The X-ray luminosity in the 0.5 — 7 keV is 
3.4 x 10 44 erg/s. The bolometric luminosity is estimated 
to be 5.1 x 10 44 erg/s. 

3.2. Radial temperature distribution 

In order to investigate the radial temperature profile, we 
accumulated spectra from four annular regions centered on 
the Gl peak. The radius ranges were chosen so that each 
spectrum contains more than 400 photons. We fitted the 
spectra with the MEKAL model to determine the tem- 



peratures, where the neutral absorption was fixed at the 
Galactic value and the metallicity of the gas was allowed to 
vary. In Figure 3(a), we plot the radial temperature pro- 
file and la error bars. The temperature in each radial bin 
was constrained with 13 ~ 20% accuracy. We found that 
there is no meaningful temperature variation with radius. 

We further restrict the spectral regions to 8" and 4" (= 
51 and 26 /i^kpc) from the Gl peak to constrain the X- 
ray emission from Gl. From the MEKAL model fitting 
to the spectra, we obtained the X-ray temperature to be 
3.7 (2.7-5.5) keV and 3.4 (2.6-4.6) keV for the regions of 
r < 4" and 8", respectively. Thus the temperature is still 
higher than 2.6 keV (90% confidence) at the center. 

From the above analysis, we found that the intracluster 
gas is consistent with being isothermal out to ~ 600 
kpc. Note that the fraction of the maximum radius to 
the virial radius, r2oo (See §6.1 for definition) is 0.41. In 
Figure 3(b), we show the 68% confidence contours for the 
temperature and the metallicity for four annular regions 
obtained from the MEKAL model fitting. We found that 
the best-fit metallicity at the cluster center is as high as 
1 solar, however, due to the large statistical errors the 
abundance gradient is not significant. 

4. IMAGE ANALYSIS 

4.1. X-ray surface brightness and galaxy distribution 

Though the original ACIS CCD has a pixel size of 
0".492, we rebinned the image increasing the bin dimen- 
sions by a factor of 4. Thus 1 image pixel is 1".968 which 
is 12.6 kpc at the cluster's frame. We restrict the energy 
range to 0.5 - 5 keV in the image analysis in order to 
avoid the hard energy band where the background domi- 
nates the total spectrum. We find that there is a second 
X-ray peak at 00:26:35.1, +17:09:38.0 (J2000). We refer 
to the second X-ray peak as G2 hereafter. The G2 peak is 
about 100 kpc southwest of the Gl peak. We investigated 
the correlation between the X-ray surface brightness and 
the member galaxies cataloged by Czoske et al. (2001). As 
shown in Figure 5a, by superposing the galaxies with red- 
shift range of 0.38 — 0.41 and V magnitude smaller than 22 
on the Chandra X-ray image, we recognized that the three 
central bright elliptical galaxies are located at positions 
that coincide with the Gl and the G2 X-ray peaks. The 
object numbers and the redshifts cataloged by Czoske et 
al. (2002) are #380 (z = 0.3936), #374 (z = 0.3871), and 
#362 (z = 0.3968) from east to west in Figure 5a. Note 
that Gl contains #380 and #374 however the X-ray peak 
position is more consistent with #380 (Figure 5a). 

4.2. 1-D fitting of the X-ray surface brightness 
distribution 

We investigate the X-ray surface brightness distribution 
by 1-dimensional and 2-dimensional model fitting in this 
and the next subsection, respectively. The latter is more 
precise in the sense that the central position of the X-ray 
emission can be included as a model parameter. Here we 
derive an azimuthally averaged surface brightness distri- 
bution centered at the strongest peak, Gl, from the 0.5-5 
keV image and evaluate the overall shape by fitting a [3- 
model, S{r) = S {1 + (r/r c ) 2 )-^ +1 /2 (Cavaliere & Fusco- 
Femiano 1976). The telescope vignetting was corrected by 
dividing the image with the exposure map. The G2 peak 
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was excluded with a circular region of 4" in radius, which 
covers more than 95% of the G2 emission, when calculat- 
ing the profile. 

We found that a single-component /3-model is rejected 
at the 99.1% level and there is clearly strong excess emis- 
sion over the model in the innermost region (Figure 4a). 
We also notice that there is a point of inflection in the 
profile at around r <~ 100 kpc, where systematic variation 
of the residuals are clearly visible. This leads us to test 
the hypothesis that the surface brightness consists of two 
components. If we introduce a double-/? model consist- 
ing of two /3-profiles with different core radii, the fit was 
significantly improved from the single /3 model and gave 
X 2 /dof = 183.2/194. In Figure 4b, we show the results 
for the double /3-modcl fitting. The derived parameters 
are listed in Table 2. Note that the background constant 
was included as a free parameter in the fit and determined 
to be C — 8.5^0 2 x 10~ 9 counts/sec/cm 2 . A strong degen- 
eracy exists between parameters r c ^ n and /?,„ (i.e. the core 
radius and the slope for the inner component of the double 
/3-model) that prevents the fit from converging properly. 
We tried different fixed values of /%„ ranging between 0.5 
and 1.5 and found that all fitting parameters except for 
r Cj i n are consistent with the result for /3j„ = 1 within the 
90% error bars. As long as 0.6 < /%„ < 1.5 the resultant 
r c ,in value is also statistically consistent with the result for 
Pin = 1. Thus we chose to fix f3 in to 1. 

We also tested a model for gas distribution in the case 
that the gas is in hydrostatic equilibrium in the NFW po- 
tential (Navarro et al. 1996). Suto et al. (1998) provided 
a useful fitting formula for the X-ray surface brightness 
distribution in the generalized form of the NFW-typc po- 
tential, p(x) oc l/(x a (l + x) 3 ~ a ), which gives a good fit for 
1.0 < a < 1.6. We refer to the formula (i.e. equation 29- 
32 in Suto et al. 1998) as the NFW-SSM model hereafter. 
It is also known that the gas density profile for the NFW 
model is well approximated by the /3-modcl (Makino et al. 
1998), though it has a slightly steeper slope in the inner- 
most region. Thus it is worth quantifying the parameter 
a. We first fitted the radial surface brightness distribution 
with the NFW-SSM model for a = 1 (i.e. the NFW case). 
The result is presented in Table 2. We found that it is 
not accepted at the 90% confidence level. If the a value is 
allowed to vary within the range 1.0 to 1.6, we found that 
the a = 1.6 case resulted in the minimum % 2 value of 199.8 
for 196 degrees of freedom. Thus the central X-ray profile 
appears to be much steeper than that expected from the 
original NFW profile. However the model with a = 1.6 
is unlikely because r s became unrealistically large (more 
than one order of magnitude larger than the value for the 
a = 1 case). 

Thus from the 1-dimensional analysis, the X-ray spa- 
tial distribution is not sufficiently described either by the 
single-component /?-model or by the NFW-SSM model and 
is significantly better fitted with the double-/? model. 

4.3. 2-D fitting of the X-ray surface brightness 
distribution 

In order to determine the X-ray emission profile of the 
ICM more precisely, we fitted the 2-dimensional surface 
brightness distribution with a model consisting of three 
(3 profiles which we consider to represent emission from 



two compact components associated with the Gl and the 
G2 peaks and "cluster" component and the constant back- 
ground; 

S(r) = J2 So,i(l + (r/r c ,) 2 )- 3 ' 3 * +1/2 + C. (1) 

We fitted the image of the central (100 image pixel) 2 = 
(3'.3) 2 region using the maximum- likelihood method. The 
center positions of the first two compact components were 
fixed at the Gl and the G2 peaks, respectively, while for 
the third component, which we consider describes the ICM 
emission, the position was allowed to vary. Because the 
slope parameters for the Gl and the G2 components were 
insensitive to the fit, we first assumed King profiles (King 
1962), namely f3 models with /3i = /?2 = 1- The core ra- 
dius of G2 was not resolved with the spatial resolution of 
the current image analysis, thus fixed at r 2 = 10 kpc as 
is typical for an elliptical galaxy. The background level, 
C, was fixed at the value that was estimated from the 
1-D analysis, 8.5 x 10~ 9 counts/sec/cm 2 . We performed 
the fit with the SHERPA package, where we included the 
exposure map to convolve the model image with the tele- 
scope and detector responses. The exposure map was cal- 
culated at an energy of 0.8 keV, representative of the spec- 
tral energy distribution. The results of the fits are shown 
in Table 3 and Figure 5. In order to check the goodness 
of the fit, we rebinned the image into two single dimen- 
sional profiles in the right ascension and declination di- 
rections (Figure 5d) and calculated the \ 2 values between 
the model and data profiles to find they are sufficiently 
small (x 2 /dof = 112.4/99 and 108.2/99 for the a- and 
(5-dircctions, respectively). The best-fit cluster center po- 
sition is significantly shifted from the Gl peak by 12". 5, 
that is 80 h§ kpc, towards the southwest direction. There 
is no optical counterpart present at this position. 

In Table 3, the results of 2-dimensional image fitting are 
listed. For Gl we obtained a core radius of 52+ g kpc. The 
cluster component was found to be much more extended 
than the Gl component, being described by a spherical 
/3-model with r c , 3 = 210^ kpc and /3 3 = OJltom- The 
90% confidence contour of the fa and r c .3 are shown in 
Figure 6. These two values are consistent with those we 
obtained from the 1-D fitting even though the cluster cen- 
ter is significantly displaced from that assumed in the 1-D 
analysis. We consider that this is because the displace- 
ment is small in comparison to the cluster core radius, 
inside which the surface brightness is flat. The small core 
radius of the Gl emission is consistent with the result from 
the ROSAT HRI though it gave a smaller [3 value. In fact 
we obtained a similarly small (3 in the 1-D analysis (see 
§4.2). This may be naturally explained by the existence 
of emission from the larger core component. In addition, 
though we fixed the j3\ value at 1 for Gl because the in- 
ner slope parameter is not sensitive to the fit, as for the 
1-D fitting, we confirmed that the resultant /3-model pa- 
rameters for all three components do not change within 
statistical errors if 0.6 < (3\ < 1.5. 

We estimated the luminosities for the three components 
by dividing the total luminosity obtained from the spec- 
tral fit according to the ratio of photon counts between the 
three. Here we assumed isothermality for the three com- 
ponents. The results are also shown in Table 3. The lumi- 
nosity of the "cluster" component is 4.5 x 10 44 erg/s and 
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dominates the total luminosity. We consider that the emis- 
sion of G2 can be attributed to an elliptical galaxy because 
our derived luminosity is within the scatter of other ellip- 
ticals (Fabbiano et al. 1992; Matsushita 2001) and there is 
a good positional coincidence with the member elliptical 
galaxy, #362. On the other hand the luminosity of the 
Gl component is higher for a typical elliptical galaxy and 
the core radius is also comparable to that of other cluster 
with a small core radius (Ota & Mitsuda 2002). We will 
discuss the properties of the Gl component in more detail 
later. 

Furthermore, we tested the significance of the ellipticity 
of the cluster image by substituting the third component 
in equation 1 for the elliptical /3-model. We found the el- 
lipticity is not significant; the 90% upper limit was e < 0.2. 
On the other hand, Soucail et al. (2000) and Bohringer et 
al. (2000) reported that the HRI image is slightly elon- 
gated in a northeast-southwest direction and can be fitted 
with a ellipse with an ellipticity of ~ 0.2 — 0.3. Soucail et 
al. (2000) also mentioned that there is a significant twist 
of the position angle at about 100 kpc from the center. 
We consider that the significant shift of the center posi- 
tion of cluster emission by 80 kpc and the presence of the 
second peak at 100 kpc off of the strongest peak in the 
Chandra image can account for the moderate ellipticity 
and the twist derived by the elliptical model analysis with 
the HRI. 

5. MASS ESTIMATION AND COMPARISON 

In this section we will derive the projected cluster mass 
profiles based on the results of the X-ray analysis and es- 
timate the cluster mass enclosed within the arc radius. 
We will directly compare these with the results from the 
lensing mass reconstruction. 

5.1. X-ray mass for the f3-model 

From the spectral and spatial analysis with Chandra, we 
found that the cluster gas is consistent with being isother- 
mal and can be described with the spherical /3-model after 
removing the local emission of Gl and G2. In what fol- 
lows, we regard the "cluster" component, i = 3 in the 
three /3-model fit, as diffuse emission coming from the 
ICM confined in the cluster dark matter potential. Since 
the Gl core radius is much smaller than the arc radius 
i r c,i *C r arc ), while the cluster core radius r c .3 ~ r arc , the 
Gl component is thought to be a minor contributor to 
the X-ray mass estimation. The effect of including the Gl 
mass will be discussed later in this section. 

Assuming the gas is isothermal and spherically dis- 
tributed and in hydrostatic equilibrium, the density of 
matter at a radius r is estimated from the /3-model as 



profile, which can be written as follows (Ota et al. 1998). 
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The projected X-ray mass density profile is given by 
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In Figure 7 we show the mass density profile calculated 
for /3-model. The cylindrical cluster mass within a certain 
radius R is calculated by integrating the projected density 
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Then from the results of 2D fitting for the cluster compo- 
nent we estimated the X-ray mass within the arc radius of 
i? arc = 35" = 220 kpc to be 

M Xl/3 (220 ft^kpc) = 0.84t°; 2 ° x 10 14 ^Mq. (5) 

Because the statistical errors of /3 and r c are coupled with 
one another (Figure 6) though kT is independently deter- 
mined from the /3-model parameters, we determined the 
error associated with the X-ray mass by evaluating it in 
the statistically allowed domain of the three dimensional 
parameter space. 

On the other hand, adopting the arc redshift of z arc = 
1.675, the strong lensing mass is estimated to be 

M lcns (< 214 h^kpc) = (3.117±0.004) x 10 14 h^M Q (6) 

from equation (2) in Tyson et al. (1998). Broadhurst 
et al. (2000) derived the strong lensing mass as Mi ons (< 
200 h^kpc) = (2.22±0.06)xl0 14 /i 50 1 M Q , which is smaller 
by about 30% than equation 6. However, comparing Mx,p 
to Mi ens , a mass discrepancy of a factor of 3 is evident. 

As shown in Figure 7, since the surface mass den- 
sity of the cluster component obtained by the isother- 
mal /3-model is less than the critical surface mass density, 
Sent = {c 2 /AnG)/{D s /D d D ds ) = 2.1 x 10 3 M Q /pc 2 every- 
where, the existence of the gravitationally multiple images 
in this cluster means that the mass distribution of the clus- 
ter is far from the isothermal /3-model and/or the existence 
of an extra mass component along the line of sight toward 
the cluster central region is required. 



5.2. X-ray mass for the NFW-model 

The X-ray surface brightness profile of the NFW po- 
tential (Navarro et al. 1996) is similar to that of the [3- 
model and can be converted from the /3-model parame- 
ters through the relations of r s — r c /0.22 and B = 15/3 
(Makino et al. 1998). We thus derived the NFW density 
profile based on the results of /3-model fitting, 



PnfwM 



p s r s 



(7) 



where p s — kTB / AitGmr 2 . We then calculated the pro- 
jected mass density profile in the same manner as equa- 
tion 3. The mass density profile is much steeper at the 
innermost region than that of the /3-model and the cen- 
tral surface mass density can be higher than the critical 
surface mass density, however, those two models arc not 
distinguishable at the arc radius (Figure 7). We estimated 
the X-ray mass to be 

M x ,nfw(220 /i^kpc) = 0.74±g;if x 10 14 h^M Q . (8) 

This shows that the NFW model does not improve the 
enclosed mass within the arc radius. Therefore, adopting 
the NFW model can not be the main solution for the mass 
discrepancy problem of this cluster. 
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5.3. Estimation of mass associated with Gl 

We will estimate the projected mass associated with the 
Gl component under some simple assumptions for the pur- 
pose of constraining its possible contribution to the X-ray 
mass estimate. Firstly, we consider the mass of the two 
bright elliptical galaxies inside Gl. Their velocity disper- 
sions and effective radii were measured to be a = 317 
km/s and r e = 5.9 /i^~ kpc for #380 and a = 275 km/s, 
and r e = 3.4 h§ kpc for #374, respectively (van Dokkum 
et al. 1996). Then calculating the projected lensing 
mass under the Singular Isothermal Sphere model, we ob- 
tain M S is ~ na 2 /(2Gr e ) = 4.4 x 10 11 M for #380 and 
M S is = 1-9 x 10 11 Mq for #374. Thus the sum of the two 
galaxies can increase Mx by only 1%. 

Second, we will treat the Gl component as the cluster 
scale substructure embedded in or projected on the clus- 
ter. One reason in support of the idea is that the observed 
X-ray properties such as the temperature, the luminosity 
and the core radius all exceed those observed for a typi- 
cal elliptical galaxy and they arc rather close to a single 
cluster. However there is no unique way to evaluate this 
component in the Mx estimation, as it is not clear how this 
is physically related to the surrounding ICM. Thus our es- 
timate of the Gl mass supposes that it is locally in hydro- 
static equilibrium and dominates the mass out to a certain 
cutoff radius (r cut ~ lOO/i^kpc) at which the inflection 
in the X-ray surface brightness is found. Thereby we ob- 
tain M G i(100 h^kpc) = 0.51±o i" x 10 14 h^M utilizing 
equation 4 and the measured temperature within a radius 
of 8" circle centered on Gl. Adding Mgi to equation 5, 
we obtain M x ,p < 1.7 x 10 14 h^M Q . We consider this 
will give a secure upper limit of the X-ray mass estimation 
under the current isothermal /3-model analysis. The total 
gas mass associated with the Gl clump within a sphere of 
r = r cut is M gas , G1 (100 h^kpc) = 5.711 j x 10 n ^ 5/2 M o 
and has a negligible contribution to the total gravitating 
mass. Thus this upper limit is still significantly smaller 
than Mi ens by a factor of ~ 2. If we compare to the Broad- 
hurst et al.'s mass, the discrepancy reduces to about 30%. 

We drew critical lines of the lens model for a source at 
z = 1.675 where the cluster component and the Gl com- 
ponent are included (Figure 8). The surface mass distri- 
bution of the Gl component is described by the equations 
(2) and (3) with above parameters for r < r cut where r 
is the projected radius from the Gl center. The surface 
mass density beyond r cut is assumed to be zero. Other 
parameters are the same as in §5.1. Inclusion of the Gl 
component turns the isothermal /3-model into the super- 
critical and a tangential critical line appears. Thus Gl has 
a significant contribution to the lensing effect. Moreover 
we notice that r c ,i is fairly close to the core size of the 
dark matter distribution found in Tyson et al. (1998). We 
then consider that this also assures the significance of the 
existence of the Gl potential. If the lens model is able to 
explain the observed multiple images, the image appearing 
on the same side of the central image, which is in the south- 
east side from the cluster center (i.e. image B in Figure 
8), should be enclosed by the tangential critical line. How- 
ever, the inclusion of the Gl component described above 
is not enough to explain the observed lensed multiple im- 
ages because the tangential critical line can not reach the 
south-east image. 



5.4. Cluster center position 

Next we directly compare the cluster center position de- 
rived by the X-ray image analysis to that of dark matter 
profile modeled by Tyson et al. (1998). We found that 
their center coordinates, (00:23:56.6, +16:53:15) (1950) are 
shifted if comparing their Figure 2 to the latest optical co- 
ordinates of central bright ellipticals from Czoske et al. 
(2001). Thus utilizing the coordinates of galaxy #380 in 
Czoske et al. (2001) and the relative distance from #380 
galaxy to the DM center, 1".5 to the east and 4". 5 to 
the south, the DM center coordinates correspond to about 
(00:26:35.8, +17:09:39.4) in J2000. Thus we found that the 
X-ray center position (00:26:35.6, +17:09:35.2) is shifted 
to the south-southwest by 4" from the dark matter center. 
If we take into account the measurement uncertainties of 
~ 2", we consider that the shift is marginal. Moreover 
Tyson et al. (1998) derived the maximum point of diffuse 
intracluster light that is not associated with visible galax- 
ies and showed it is displaced 3" west-southwest from the 
DM center. The X-ray center position is closer to this than 
the DM center. 

6. DISCUSSION 

6.1. Properties of the intracluster gas 

We will discuss the X-ray properties of the gas focusing 
on the X-ray luminosity-temperature (Lx — T) relation, 
gas-mass fraction and the metallicity based on the Chan- 
dra analysis. 

The location of the cluster on the Lx — T plane will pro- 
vide important information for understanding the physical 
status of the ICM in the cluster central region. The ob- 
servational data of clusters are known to exhibit a large 
scatter around the best-fit power-law model. Fabian et al. 
(1994) noted that the Lx — T relation is roughly divided 
into two sequences whether the clusters are XD or non-XD. 
The XD cluster is the cluster in which the X-ray emission 
is highly concentrated in the central giant elliptical galaxy 
and the X-ray peak position coincides with the center of 
the central galaxy. In the non-XD cluster, the emission is 
diffuse and the X-ray peak position is largely offset from 
the central giant galaxy. Arnaud & Evrard (1999) derived 
the Lx — T relation with the analysis sample restricted to 
non-XD clusters and showed a smaller scatter. A different 
approach was taken by Ota & Mitsuda (2002), concerning 
the cluster core radius. They performed a systematic anal- 
ysis of 79 distant clusters with the ROSATRRl and ASCA 
to study the X-ray structure of clusters in 0.1 < z < 1. 
They determined the average X-ray temperatures and the 
bolometric luminosities with ASCA and the X-ray surface 
brightness distributions with the ROSAT HRI by utiliz- 
ing the isothermal /3-model, and found there is not any 
significant rcdshift dependence in the X-ray parameters 
including the temperature, /3-model parameters, and the 
central electron density. They discovered that the distri- 
bution of the core radius shows distinct two peaks at 60 
kpc and 220 kpc. If dividing the cluster samples into two 
subgroups corresponding to the two peaks in the core ra- 
dius distribution, they show differences in the X-ray and 
optical morphologies and in the Lx — T relation. In par- 
ticular, the normalization factor of the Lx — T relation 
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significantly differs according to the core sizes: at a cer- 
tain temperature, the luminosity is higher for a cluster 
with small core radius and approximately Lx follows r" 1 . 
From these observational results, they suggested that the 
clusters are divided into at least two subgroups according 
to the core radius. In §4, wc showed that the X-ray spatial 
distribution of CL0024+17 is described with a superposi- 
tion of two [5- model components with r C) i ~ 50 kpc and 
r C 3 ~ 210 kpc. Therefore, the main cluster component is 
classified as a non-XD cluster or a large core radius clus- 
ter, while the Gl component is classified as an XD or a 
small core radius cluster. 

We will compare our result for CL0024+17 with Arnaud 
& Evrard (1999) and Ota & Mitsuda (2002). First, em- 
ploying the Lx — T relation of Arnaud & Evrard (1999), 
the luminosity is expected to be logL^ = 44.69 ± 0.05 
for the observed temperature of 4.5 keV. Thus the X- 
ray luminosity of CL0024+17 determined with Chandra, 
logLx.boi = 44.7, is in a good agreement. Moreover, we 
compare our result to the Lx — T relation for clusters with 
large(r c > 135 kpc) core, derived by Ota & Mitsuda (2002) 
since the main cluster component whose core is r Cj3 ~ 210 
kpc is responsible for the cluster luminosity (See Table 3). 
We found that the luminosity of CL0024+17 is smaller 
than the other distant cluster samples with similar tem- 
peratures, however, it is within a scatter of the data points 
and is quite consistent with the mean relation of the large 
core (r c > 135 kpc) clusters. Therefore in terms of the 
Lx — T relation, the X-ray emission of CL0024+17 agrees 
with other clusters that do not have strong central X-ray 
emission dominating the total luminosity. Adopting the 
estimated Gl luminosity of 5.5 x 10 43 erg/s (Table 3), the 
temperature of 1.2 keV is obtained from the Lx — T re- 
lation for the small core radius clusters. The obtained 
temperature is significantly lower than that constrained 
by the X-ray observation (§3.2). This might be the sig- 
nature that the Gl component is now in the process of 
merging with a cluster or has undergone a recent merger 
and therefore the characteristics of the Gl component can 
not be estimated by using a relation for relaxed clusters. 

Since the observed X-ray spectrum is well fitted with the 
single-temperature MEKAL model (§3), the contribution 
of additional cool emission is suggested to be not impor- 
tant in this cluster. We can then estimate the time scale 
of radiative cooling at the cluster center. By deprojecting 
the central surface brightness of the (3 profile, we obtained 
the central electron density for the cluster component as 
^e0,3 = (5.8 ± 0.6) x 10~ 3 cm" 3 . Then the cooling time 
scale is 11.4 Gyr at the center. This is longer than the 
age of the Universe, 7.9 Gyr. Thus radiative cooling is not 
likely to be effective. On the other hand, the electron den- 
sity is higher for Gl, where n e o,i = (2.5±0.5) x 10~ 2 cm" 3 . 
This yields a much shorter cooling time of 2.8 Gyr. Thus 
cooling may occur very effectively in the Gl region. How- 
ever it is contradictory that we did not find a significant 
temperature decrement in the spectral analysis (See §3.2). 
A similar situation has been found in some nearby clusters 
with short cooling timescales (e.g. Tamura et al. 2001) and 
other possibilities to prevent cooling have been proposed 
(e.g. Bohringer et al. 2002). We will discuss the possibility 
of a two-phase state of the ICM in the next subsection in 
detail. 



We estimate the gas-mass fraction in the cluster within 
the virial radius in order to compare it with the mean 
baryon fraction in the Universe determined by the Wilkin- 
son Microwave Anisotropy Probe (WMAP) experiments 
(Spcrgel et al. 2003). In this analysis, we use the standard 
set of the cosmological parameters determined by WMAP, 
Q m = 0.3, Cl\ = 0.7, and h = 0.7. The virial radius is 
defined as the radius at which the cluster density is equal 
to 200 times the critical density of the universe at the clus- 
ter rcdshift, r 2 oo = 1-39 h^ 7 Mpc. Then the virial mass 
(the spherical hydrostatic mass), gas mass and gas-mass 
fraction calculated at r 2 oo under the /3-model are 

M 2O o=4.6t^xlO 14 ^iM , (9) 
M gas = 6.5t^ x 10 13 h 5 7 /2 M Q , (10) 

/gas = M gas /M 200 = 0.14±g:8S h- 3 7 /2 . (11) 

The gas-mass fraction at the virial radius is consistent 
with the Universal baryon fraction derived by WMAP, 
0;,/O m = 0.16, within 35% accuracy. Thus we consider 
the virial mass to be properly evaluated with M200 an d 
the X-ray temperature to represent the virial temperature 
of the cluster potential. 

As shown in the spectral analysis, the average iron abun- 
dance is about twice as high as the typical value. Since the 
mctallicity is a quantity defined relatively to the amount 
of hydrogen, there are two situations that will explain this 
high value: a medium rich in iron or a medium poor in 
hydrogen. Because the gas-mass fraction (Equation 11) 
is comparable to that observed in other clusters (Mohr 
et al. 1999; Ota & Mitsuda 2001), it is plausible that 
the ICM was highly enriched in iron. As for nearby clus- 
ters, the iron mass in ICM, Mp c , is known to have a clear 
correlation with the optical luminosity of E+S0 members 
(Arnaud et al. 1992) and with the total blue luminosity 
(Rcnzini 1997). Schneider et al. (1986) mentioned that 
the luminosity of CL0024+17 places it among the richest 
of any known clusters. Because the previous measurement 
of the cluster metallicity showed no significant rcdshift evo- 
lution (e.g. Mushotzky & Loewcnstcin 1997), and as the 
contribution of the metal production by type la SNe after 
z = 0.4 can be negligible according to the model calcula- 
tion by Mihara & Takahara (1994), we will compare the 
iron mass observed at z = 0.395 to the optical luminosity 
corrected for passive evolution in ellipticals below. Utiliz- 
ing the local relation, M Fo ~ l-6± 3 ;'i x IO^L^e+so ( Ar " 
naud et al. 1992) and the visual luminosity of ellipticals 
evolved to the present day Ly(z = 0) = 1.2 x 10 12 /i 50 2 L Q at 
r = 400 h^Q kpc (Smail et al. 1997), the iron mass is calcu- 
lated as 1.9+ 4 'Jx 10 10 /i5 2 ' 5 M Q , where we adopted the best- 
fit slope of the M — Lv,e+so relation, i.e. 1, and included 
only the la error of the normalization factor. Though the 
relation of Arnaud et al. (1992) was derived at a radius of 
3 Mpc, we assumed here that iron distribution is pro- 
portional to that of the galaxy and then it is applicable for 
a smaller radius. On the other hand, the Chandra analy- 
sis yields M Fo (< 400 /i^ 1 kpc) = (2.5 ±0.2) x lO lo /i^ o 2 - 5 M 
where we adopted the mean metallicity measured within 
600 h^ kpc, 0.76 solar. Because this is by about 30% 
higher than the above calculation, there may be an abun- 
dance variation from cluster to cluster, as suggested from 
the recent Chandra results of metallicity maps: the inho- 



8 



Ota et al. 



mogeneous metal distributions were found in the core re- 
gions of some nearby clusters such as Perseus (Schmidt, 
Fabian, & Sanders 2002) and A2199 (Johnstone et al. 
2002). However, since the iron mass of CL0024+17 is 
consistent with the Mp c — iy.E+so relation within the 
lcr errors, the variation is not statistically significant. To 
confirm this, we need future observations of the metallicity 
map with higher sensitivities. 

The observed high metallicity can be explained by metal 
ejection from the elliptical galaxies. Since the spiral galaxy 
fraction, / sp = 0.40, is not different from other clusters 
with a similar mass (Smail et al. 1997), the above result 
indicates that very effective galaxy formation had occurred 
in the central region of CL0024+17. Besides, the optical 
observations showed that the spatial distribution of galax- 
ies is centrally concentrated and has a very compact core 
(Smail et al. 1996; Schneider et al. 1986). Then because 
of the high concentration of bright elliptical galaxies the 
metal will also follow the concentrated distribution. As 
mentioned above, there may also be a complex metal dis- 
tribution in the cluster core. However due to the poor pho- 
ton statistics, we are not able to derive the 2-dimensional 
abundance map. For this reason we studied the radially- 
averaged abundance profile in §3.2. Our result shown in 
Figure 3b indicates that the metal abundance is as high as 
~ 1 solar at the innermost region and then not in conflict 
with an abundance profile in which Zp c ~ 1 within the 
central r < 10" region and Zp e ~ 0.3 at the outer region. 
In order to further constrain the spatial distribution and 
the history of the metal production, better photon statis- 
tics for the X-ray spectrum are required. 

6.2. Mass discrepancy problem 

We found from the Chandra observation that the X-ray 
surface brightness distribution is mainly described by a su- 
perposition of two extended components well fitted with 
spherical /3-profiles. Furthermore we do not find any sig- 
nificant temperature structure. Thus it is likely that the 
gas is relaxed in the cluster potential and, therefore, that 
hydrostatic equilibrium will be a good approximation in 
the X-ray mass estimation. We further discuss any possi- 
ble cause of the mass discrepancy below. 

As shown in the previous section, we can consider that 
the X-ray temperature represents the virial temperature of 
the cluster component since the gas mass fraction of the 
cluster component within the virial radius is consistent 
with the universal baryon fraction obtained by WMAP. 
Therefore, the mass discrepancy reported in the previous 
section is telling us that there is a lack of our current un- 
derstanding for the nature of the cluster central region. 

In some nearby clusters, particularly XD clusters (For- 
man & Jones 1990), an additional cool component is re- 
quired to explain the X-ray emission of the central ~ 
100 /i^kpc regions (Makishima et al. 2001). On the con- 
trary Makishima et al. (2001) also suggested that non-XD 
clusters appear to be isothermal with little metallicity gra- 
dient toward the center. Taking into account the fact that 
the current target cluster is a non-XD cluster and that the 
total X-ray emission is not dominated by the central Gl 
emission, our results shown in §3 seem to be more consis- 
tent with their picture of non-XD systems. However if the 
ICM of CL0024+17 is in a two-phase state and the thermal 



pressure of the gas in each phase balances against gravita- 
tion, it would greatly increase the X-ray mass estimation. 
Thus we attempted to fit the spectra for both (1) the clus- 
ter region outside 100 kpc from Gl, i.e. 0'.26 < r < 1'.5 
and (2) the Gl region, r < 0'.26 with a two-component 
MEKAL model. Because we are not able to constrain the 
two temperatures simultaneously under the current pho- 
ton statistics, we fixed the temperature and the metallicity 
for the cool phase respectively at 1 keV and 1 solar while 
the spectral normalization was adjustable. For the region 
(1), we found that the flux of the cool phase is only about 
4% of the hot phase and the temperature of the hot-phase 
is 5.5 (4.4 - 7.2) keV, which is consistent with the result 
of the single-phase model shown in Table 1 within errors. 
Thus the spectrum is well represented by the single-phase 
model. This is also consistent with the results of nearby 
non-XD clusters. For the region (2), because the radiative 
cooling time scale of the Gl emission is shorter than the 
age of the Universe (See §6.1), the gas might be in a two- 
phase state. However, from the two-phase model fitting, 
we found that the upper limit of the cool-phase flux is less 
than 1% of the hot-phase and the hot-phase temperature, 
4.2 (3.5 - 5.3) keV, is again consistent with that obtained 
for the single-phase model. Thus we consider that it is 
unlikely that the Gl gas is in a two-phase state and that 
the virial temperature is much higher than that derived 
from the single-phase model. 

Czoske et al. (2001, 2002) measured the redshift distri- 
bution of galaxies in the direction of CL0024+17 and re- 
vealed the presence of foreground and background groups 
of galaxies that align along the line of sight. Czoske et al. 
(2002) suggested from their results that there was a high 
speed (~ 3000 km/s) collision between the CL0024+17 
cluster and a second cluster with a mass about half that 
of the main cluster ~ 3 Gyr ago. They suggested that tak- 
ing into account the projection effect of the second cluster 
mass may solve the mass discrepancy problem. They also 
mentioned that during the collision the X-ray gas would be 
highly disturbed due to propagation of shock waves (e.g. 
Takizawa 1999) and hydrostatic equilibrium of the gas 
component will break down, however, after several Gyrs 
the gas will settle down to an equilibrium state. Roet- 
tingcr et al. (1996) noted that, based on numerical simu- 
lations, the temperature structure is one of the strongest 
indicator of recent merger activity, however, we do not 
find any substantial temperature structure in the observa- 
tional data as shown in Figure 3(a). We also found from 
the spectral analysis of the Gl peak that the temperature 
is still higher than 2.6 keV at the central r < 4" region. 
Thus our results show that the cluster is in the state of 
several Gyr after the merging, and the gas should have 
settled down to an equilibrium state and then trace the 
underlying dark matter potential, which is supportive of 
Czoske et al. (2002) 's scenario. 

As discussed above, the total cluster mass is considered 
to be properly evaluated from the X-ray observations at 
the virial radius based on the consistency of the gas-mass 
fraction with the Universal baryon fraction. On the other 
hand, the large mass discrepancy detected in the central 
region may indicate that there was a merging along the 
line of sight which disturbs the mass distribution in the 
cluster core region, as suggested by Czoske et al. (2002). 
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The Gl peak that we found in the X-ray image may be 
related with the merging. In addition if we are seeing the 
two collided cores superposed along the line of sight, the 
mean rcdshift values may be slightly different between the 
cores. However we could not constrain the redshifts of the 
Gl and the surrounding region from the X-ray spectral 
data under the current limited photon statistics and the 
energy resolution for the Fe line emission. We expect that 
the future ASTRO-E2 mission will constrain the line-of- 
sight structure of the cluster mass distribution with its 
finest spectral resolution. 

Although the previous precise lens modelings (Tyson et 
al. 1998; Broadhurst et al. 2000) were performed with an 
a priori assumption that CL0024+17 is an isolated single 
cluster of galaxies, our results have provided new evidences 
which show that this is not the case. One of the strongest 
pieces of evidence is the consistency of the gas-mass frac- 
tion of this cluster obtained from our X-ray results with 
the universal baryon fraction. In the case that the previ- 
ous lens models (Tyson et al. 1998; Broadhurst et al. 2000) 
describe the cluster mass distribution up to the virial ra- 
dius correctly, the gas mass fraction falls to 1/3 of the 
universal baryon fraction. Thereby it will be a key sub- 
ject to construct a new lens model based on the current 
X-ray results and the results of the optical redshift survey 
by Czoske et al. (2002) which suggests the existence of the 
merger in the line of sight, and to examine the details of 
the physical nature of the Gl component by further X-ray 
observations. 

7. SUMMARY 

With Chandra we have performed the spectral and the 
spatial analysis of the lensing cluster CL0024+17. The 
temperature profile is consistent with being isothermal 
and the average temperature is ~ 4.5 keV. We detected 
a strong redshifted iron line with a corresponding iron 
abundance of ~ 0.8 solar, which is one the largest val- 
ues among known clusters. We found that the radial X- 
ray surface brightness distribution is not fitted with a 
single-component /3-modcl or an NFW-SSM model and 



that there is significant excess emission centered at central 
bright elliptical galaxy. The X-ray surface brightness dis- 
tribution was well described by a 2-dimcnsional /3-model of 
two extended components: the Gl component with a small 
(<~ 50 kpc) core radius and the surrounding main cluster 
component with a core radius ~ 210 kpc and whose cen- 
ter is shifted by ~ 80 kpc from the Gl center. We derived 
the X-ray mass within the arc radius assuming hydrostatic 
equilibrium and compared it to the strong lensing mass. 
Our result clearly showed that the X-ray mass is signifi- 
cantly smaller by a factor of 2-3. We then discussed the 
possible cause of the mass discrepancy concerning the two- 
phase spectral model, the cluster merging and the effect 
of the secondary potential in the lensing effect. Although 
we do not rule out the possibility that Gl is a remnant of 
a cluster merger, because of the absence of any substan- 
tial temperature structure indicative of a recent merger 
the gas seems to have relaxed in the cluster potential. It 
is also true that Gl plays an indispensable role to make 
the surface mass density of the lens supercritical, however 
it is not sufficient to reconcile the large mass discrepancy. 
Considering the fact that the X-ray measurement of the 
gas-mass fraction at the virial radius is consistent with the 
universal baryon fraction, we suggest that it is important 
to further study the physical nature of the cluster core and 
clarify its relation to the merging process, and incorporate 
it into constructing a new lens model together with the 
latest optical information of the galaxy distribution. 
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Fig. 1. — X-ray image of CL0024+17 in the 0.5 — 7 keV band obtained with the Chandra ACIS-S3 with overlaid contours. The image is 
adaptively smoothed and corrected with the exposure map. The background is not subtracted. The image and the contours are logarithmically 
scaled. The point sources detected with the ROSAT HRI, S1-S5 (Soucail et al. 2000) are shown with the ROSAT error box of 16" X 16". 
The boundary of the ACIS-S3 CCD chip is shown with 8' X 8' box. 
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Fig. 2. — Chandra ACIS-S3 spectrum of CL0024+17 (r < l'.5) fitted with the MEKAL model. In the upper panel, the crosses denote the 
observed spectrum and the step function shows the best-fit model function convolved with the telescope and the detector response functions. 
In the lower panel, the residuals of the fit in units of a arc shown. 
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Fig. 3. — (a) Radial temperature profile determined from the MEKAL model fitting to the four annular regions. The vertical error bars 
are la. (b) Confidence contours of the temperature and the metal abundance for the same four annular regions used in the left panel. The 
single-parameter 68% error domains are shown as contours with the thick-solid, thin-solid, thin-dashed, thin-dotted lines from the inner to 
the outer annuli. 
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Fig. 4. — Results of the 1-dimensional surface brightness distribution fitting with the single /3 model (a) and the double fi model (b). In 
each panel, the crosses show the observed surface brightness and the solid line shows the best-fit model. The constant background is shown 
with the horizontal dotted line. For the double fi model, the inner and the outer components are respectively shown with the dashed and the 
dash-dotted lines. 

Fig. 5. — Results of the 2-dimcnsional surface brightness distribution fitting with the three (3 profiles, (a) Chandra ACIS-S3 image of 
the central 3'. 3 X 3'. 3 region of CL0024+17 in the 0.5 — 5 keV energy range. The small circles denote the positions of the galaxies with 
0.38 < 2 < 0.41 (Czoske et al. 2001). A closer view of the central 30" X 30" region is also shown in the upper-left panel, (b) The best-fit 2D 
image of the three /3-models, overlaid by the contours with logarithmic spacing, (c) Residuals of the 2D image fitting, (d) Residuals projected 
to the x- and y- directions in units of a. 
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Fig. 6. — x 2 contours of the 2-dimcnsional image fitting. The 90% single-parameter error domain for /?3 and r Cj 3 of the cluster component 
is shown. The position of the \ 2 minimum is denoted with a cross. 



Fig. 7. — X-ray mass density profiles for the /3-model (the thick line) and the NFW model (the thin line) of the main cluster component. 
The 90% error ranges are shown with the dashed and the dotted lines, respectively. The arc radius of 220/i,T 1 kpc is denoted with the arrow. 
The projected dark matter density profile derived by Tyson et al. (1998) was plotted with the solid curve using the fitted parameters given 
by Shapiro & Ilicv (2000). The critical surface density, X C rit> is indicated with the horizontal dash-dot line. 
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Fig. 8. — Critical lines for a source at z = 1.675 are overlaid on the B band image of CL0024+17 obtained with CFHT (the original image 
was described in Czoske et al. (2001)). North is up and east is left. The multiple arc images are marked with A-E in the same manner as 
Colley et al. (1996). No critical line appears in the case of a cluster lens for the main cluster potential (see §5.1). On the other hand, when 
adding the Gl potential the critical lines appear (the inner and the outer curves are the radial and the tangential critical lines respectively), 
however, it is not enough to reach to the southeast image (see §5.3). The object numbers of the central elliptical galaxies in the catalog of 
Czoske et al. (2001) are also shown. 
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Table 1 

Results of the MEKAL model fitting to the overall spectrum 



Parameter Value (90% error) 

iV H (cm" 2 ) 4.2 x 10 2U (F) 

kT (keV) 4.47 (3.93 - 5.30) 

Abundance (Z Q ) 0.76 (0.45 -1.13) 

Redshift 0.395 (F) 

Normalization* 5.72 (5.20 - 6.25) 

X 2 /dof 47.5/45 
/x,o.5-7 (erg/s/cm 2 ) b 4.1 x lO" 13 

L X;0 .5-7 (erg/s) c 3.4 x 10 44 

L x ,boi (erg/s) d 5.1 x 10 44 



"Normalization factor for the MEKAL model, J n e nudV / Att D\(l + z) 2 in 10 18 cm 5 , where Da is angular size distance to the cluster. 
ft X-ray flux in the 0.5-7 keV band. 
c X-ray luminosity in the 0.5—7 keV band. 
d Bolomctric luminosity. 
(F) Fixed parameters. 



Table 2 

Results of the 1-D image fitting 



Model 


(component) 


cts/s/cm 2 /(h 50 1 kpc) 2 
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"IKo kP c 


X'fdof 


Single-/? 
Double-^ 


(inner) 
(outer) 


1.6+^ x 10- 9 
2.8 +0 8 7 x 10~ 9 
9.6±°; 2 x 10- 10 


u - oo -0.04 

1.0(F) 

66+ 08 

u - oo -0.06 


17.2lt;^/109+^ 

7.0li' 4 4 /45ti° 

29.5l 6 5 ;° /187t 3 3 « 


242.5/196 
183.2/194 


Model 
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cts/s/cm 2 /(h 50 1 kpc) 2 
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"/Ko kP c 


X'/dof 


NFW-SSM 


1(F) 


1.8+^ x 10- y 


o O+0.8 

o-J-0.6 


65±iS/415±^ 


229.3/196 



(F) Fixed parameters. 

"Central surface brightness of the /3-profile or the NFW-SSM model in the 0.5 — 5 keV. 



Table 3 

Results of the 2-D image fitting with the three /3-models 



Model 


i 


Center position 






ft 






component 




RA,Dcc in J2000 


cts/s/cm /(h 


r,n lk P c ) 2 




" l h J2 k P c 


erg/s 


Gl 
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00:26:36.0,+17:09:45.9(F) 


3-ito^xio- 


-9 


1 (F) 




5.5 x 10 43 


G2 

Cluster 
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00:26:35.1, +17:09:38.0 (F) 
00:26:35.6, +17:09:35.2 b 


4.7+J' 5 x 10" 
8-7±oi| x 10- 


-9 
-10 


1 (F) 
71+ - 07 

u - ' -0.06 


1.6/10 (F) 
32.8^4 7/2101^ 


3 x 10 42 
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(F) Fixed parameters. 
"Central surface brightness of the /3-profile in the 0.5 — 5 keV. 
6 The 90% errors are ±1".3 for RA and ±1".5 for Dec. 
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